Interaction of Vortices in Complex Vector Field and Stability of a "Vortex Molecule" 
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We consider interaction of vortices in the vector complex Ginzburg-Landau equation (CVGLE). In 
the limit of small field coupling, it is found analytically that the interaction between well-separated 
defects in two different fields is long-range, in contrast to interaction between defects in the same 
field which falls off exponentially. In a certain region of parameters of CVGLE, we find stable 
rotating bound states of two defects - a "vortex molecule" . 
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The complex Ginzburg-Landau equation (CGLE) is 
a paradigm model for qualitative description of weakly 
nonlinear oscillatory media (see for review |l]]). This 
equation is a generic form which is obtained as the ampli- 
tude equation in the vicinity of a Hopf bifurcation. There 
is a vast literature on vortex solutions of CGLE and on 
interactions and instabilities of vortices [j]J^| . 

Recently, the attention has been brought to the vector 
extension of CGLE. The problem of nonlinear dynamics 
of a 2D complex vector field arises most naturally in the 
context of nonlinear optics, where the order parameter 
is the electric field envelope in the plane normal to the 
direction of propagation |^-|| ; the fields A± can be iden- 
tified with the two circularly polarized waves of opposite 
sense. Alternatively, coupled complex fields can be inter- 
preted as amplitudes of interacting nonlinear waves Q . 
A distinguished feature of the vector Ginzburg-Landau 
equation (VGLE) is a possibility of transition between 
two "phases" , which can be characterized by either "mix- 
ing" or "separation" of two "superfluids" . Defects (vor- 
tices) can exist in both "superfluids", and transitions 
between alternative core structures are possible in 
this sense, VGLE could be viewed as a toy model of 3 He 
dressed down to two dimensions. 

A particularly intriguing possibility, suggested in 
Ref. [pj, is the formation of a bound pair of defects in 
the two fields, i.e. a vortex "molecule" with dipole struc- 
ture. We shall show in this Letter that such a "molecule" 
cannot in fact exist in the model with real coefficients, 
but is readily formed in the vector model with complex 
coefficients (CVGLE). The latter form appears, in partic- 
ular, as an amplitude equation near the lasing transition 
j|[7) . Recent simulations of CVGLE j| have shown spiral 
wave patterns with an exceptionally rich structure where 
both separated (but closely packed) defects in the two 
fields and "vector" defects with a common core (called 
argument and director vortices in Ref. Q] ) could be seen. 
There has been so far, however, no theoretical studies of 
isolated defects in CVGLE and their interactions. 

In this Letter we study interaction of defects in CV- 
GLE in the limit of small coupling between two complex 
fields. We have found that the interaction between a well- 
separated pair of defects in two different fields is always 



long-range (power-like), in contrast to the interaction be- 
tween defects in the same field which falls off exponen- 
tially as in a sing le CGLE §. In a certain region of 
parameters of CVGLE we found stable rotating bound 
states of two defects - a "vortex molecule". Analytical 
results are in excellent agreement with simulations. 

Under appropriate scaling of the physical variables, the 
vector equation for two symmetric interacting complex 
fields A± acquires the universal form 



d t A ± =A ± -(l + ic)(\A ± \ 2 + g\A T \ 2 )A, 
+ (l + ib)W 2 A ± , 



(1) 



where real parameters b and c are, respectively, the ratio 
of dispersion to diffusion and the ratio of conservative 
to dissipative nonlinearity, and the complex parameter 
g = g r + igi characterizes the magnitude of the coupling. 

A scalar defect of Eq. (|]) with unit topological charge 
is a one-armed spiral in A + (or A-) field, while a vec- 
tor defect is formed by two scalar defects with a common 
core. Simple energy considerations, applicable in the case 
of VGLE with real coefficients, point out that scalar de- 
fects in different fields tend to separate or unite, respec- 
tively, when g is positive or negative. For the CVGLE, 
no energy integral can be defined, but the limit g — > is, 
of course, distinguished, since the interaction ceases, and 
defects in both fields, described by the usual Hagan's |ic[ | 
rotating spiral solution, are mutually independent. 

Defects in two different fields. At \g\ -C 1 the interac- 
tion can be treated perturbatively. For g = the scalar 
defect of Eq. ([!]) in either field centered at the origin is 



A(r, 9) = F(r) exp i[-wt + 9 + ip(r)], 



(2) 



where (r, 9) are polar coordinates, ui = c + (b — cjkg, are 
the rotation frequency, and fco is an asymptotic wavenum- 
bcr emitted by the spiral. The real functions F and 
tjj have the following asymptotic behavior for r — > oo: 
F — > yl — fcp and x ~~ * &0> where x = V''( r ) is a 
local wavenumber. The wavenumber fco is determined 
uniquely for given 6, c [|l0| . 

A form more convenient for the analysis, obtained after 
setting A± = u±y/(l + uJb) /(l + bc)e~ lu)t and rescaling 

V -> VV(l + M/(l + & 2 ), &->(! + bu)8t, is 
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(1 — ib)dtu± = (1 + iQ)u± 

- (l + iq)(\u±\ 2 + g\u T \ 2 )u ± + V 2 u ± , (3) 

where q — (c — b)/(l + be), ft = (ui — c)/(l + wc). 

Due to the interaction, the positions of defects r±(t) 
become slow functions of time, so that the instantaneous 
drift velocity r± = v± = 0(|<?|)- Rewriting Eq. (||) in 
the comoving frame, we obtain in the first order 



(1 - ib)v± ■ Vw± + (1 + i£l)u± 

(1 + iq) (\u ± \ 2 + g\u T \ 2 ) «± + V 2 u± = 0, 



(4) 



The imaginary part of the advective term can be ab- 
sorbed by transforming u± — > u± exp (ii6v± ■ x) , which 
accounts for the Doppler shift in the emitted wave pL 
We concentrate upon one defect, say, that marked by 
the index + (which we will further omit), take its position 
as the origin and view interaction with its counterpart as 
a perturbation. Expanding in g, we write 

u± = (-F(|r - r±|) + w±) expi[9± + ^(|r - r±|)], (5) 

where r± are positions of the zeroes of the respective 
fields and 6± are polar angles about these points; w± 
is an 0(|g|) correction. Substituting the ansatz (|1) into 
Eq. (||), we obtain the first-order equation H + = 0, 
containing the linear operator (cf. Q) 



2i 

H=-(l + iq)F 2 (w + w*) + Aw + -^d e w 4 
and the inhomogeneity 
* = -(! + iq)gFF 2 (f) + (F' + iF X )w ■ n - 
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Here r = |r — R|, R = r — r + is the separation between 
the defects, n = (cos 6, sin 6) is the unit vector along R, 
and A = V 2 - F~ 1J! 7 2 F. 

The operator H has two Goldstone modes correspond- 
ing to the transnational symmetry in the plane. Therefore 
the inhomogeneous equation H + ^ = has bounded so- 
lution if ^> is orthogonal to the adjoint zero modes w + of 
the homogeneous problem. Thus, the drift velocity can 
be derived from the solvability condition (see Ref. B) 



Re J w + (r)*(r)d 2 r = 0, 



(0) 



The operator Ti. is not self-adjoint for any q ^ 0. As a re- 
sult, the adjoint zero modes cannot be expressed through 
the translational modes Vuo, and have to be computed 
directly by solving the equation Tt + (w + , w + *) = 0. The 
adjoint operator 7i + is 



TL 



-F 2 (w^ 



,+* 



iq{w — w )) + Aw 



2l a h 



2i 

—d r {r X Fw+) 



(7) 



We take note that the Goldstone modes of Eq. (||) are 
first harmonics w = e"^ 6,+ ^- > Vuo. In view of the orthogo- 
nality of the cigenfunctions of the operator and its adjoint 



with distinct eigenvalues, the zero modes of the adjoint 
operator must contain the first harmonics as well. To 
solve the equation TC + = 0, we first separate the real and 
imaginary parts of w + = t + + is + and then take the lat- 
ter 's first harmonics (t+ , s+) = (T+ (r) , S + (r)) e ie + c.c. 
Thus, we derive from Eq. ([ 
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0, 



(8) 



- F- 1 \7 2 F. A typical 



where A r = d 2 + r~ 1 d r — 
structure of the adjoint mode in shown in Fig. 1. One 
can find that the adjoint mode decays exponentially for 
r — > oo (see also Ref. Q). 




FIG. 1. The real and imaginary parts (solid/dashed lines) 
of T + and S + components of the adjoint mode for q = 1.3. 

The solvability condition (^|) can be rewritten in a com- 
pact form by defining the complex velocity v — v x — iv y , 
so that v • n = Re(i>e* e ), v x n = Im(we* e ). Assuming 
that R is directed along the x axis, v x and v y coincide 
with the radial and tangential velocity components, re- 
spectively, v r and v T . Expressing also the first harmonics 
of the cross-coupling term through the angular integral 
F(r) = tt^ 1 Jq cos 9 (F(f)) 2 d9, we derive from Eq. (||) the 
equation for the complex velocity 



v = X lf -f r 



T 



FF(r)rdr. 



(9) 



where 7 = j r + iji = (1 + iq)g, and the friction factor is 



1 = 



F'(r)T+ + X F 



iF 



S+ 



rdr. 



This expression is the principal analytical result of our 
work. One can see immediately that, with a generic com- 
plex T + (r) and S + (r), the velocity can be modified sim- 
ply by rotating the argument of the complex interaction 
parameter g. In this way, one can ensure that the radial 
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component v r = Re(v) vanishes at some finite R at least 
in a certain interval of arg(<?) . If this equilibrium position 
is stable, it signals the formation of a bound state, which 
we call a vortex molecule. Since at the equilibrium dis- 
tance the tangential component of the velocity given by 
v T = Im(zj) is generally non-zero, this bound state must 
rotate with a certain angular velocity. 

The bound states cannot form in the real VGLE. In 
this case the operator TC is self-adjoint, and the ad- 
joint mode is just the translation mode, (T + ,S + ) = 
(F' ,—iF/r). From Eq. (^|) one derives then v r ~ 
g J FF'Frdr and v T = 0. Since F(r) is a monotonic 
function, v r does not change sign. One finds that the 
vector defect is unstable for g > and stable otherwise. 
The eigenvalue A of the mode responsible for splitting the 
vector defect as function of g is shown in Fig. 4 (inset). 

For R > 1 and R < 1 Eq. (|) can be calculated 
analytically. Using the asymptotic expansions F 2 (f) ~ 
1 — k 2 + k /(qf) + ... valid for R ^> 1 we obtain 



Jo 



cos 6d6 



Vr 2 + R 2 - 2rRcos( 



kpr 
2qR 2 ' 



(10) 



For i?< 1 one uses F(f) = F(r) +F(r)'(R/r-R cos 6), 
leading to F = -RFF'(r). The solvability condition (|) 
yields then v = —R^ai for R -C 1, v = R~ 2 ~f r a2 for 
R^> 1 with the constants a% 2 given by 



Oi\ = 



a 2 



I- 1 F 2 F'(r) (t+ + ^S+^j rdr, 



^0 -r-l 
2q 



I- 1 r F(r) ( T+ + ^S+ 

JO V lr 



r 2 dr. 



(11) 



integrals in Eq. (^) and calculated both the constants 0:1,2 
(Fig. 2) and radial and tangential velocities as functions 
of the distance R (Fig. 3). We find that for = 0, g r > 
and q < q c « 0.52 the defects repel one another (v r < 0) 
at small distances, so that a vector defect is unstable, as 
in the case of real VGLE. Surprisingly, for larger q > q c 
the real part of a\ changes sign, and the defects bind at 
R = forming a stable vector defect. Correspondingly, 
for g r < the vector defect is stable at q < q c as for real 
VGLE, but becomes unstable at q > q c . 
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FIG. 3. The radial and tangential velocities v r , v T as func- 
tions of the separation distance R for q — 1.3 and g = —0.01. 
The solid line has been computed using Eq. (^|) , circles present 
results of simulations of Eq. d3J) . 




FIG. 2. The real (solid line) and imaginary (dashed line) 
parts of ai and Q2 as functions of g (for g real). 

Provided F(r),x and the adjoint mode T + ,S + are 
known, the dependence of velocity on distance R can be 
found explicitly. We used a matching-shooting algorithm 
to find a stationary spiral solution and a corresponding 
adjoint eigenmode. Then, we evaluated numerically the 




FIG. 4. Equilibrium separation R vs. q obtained from 
Eq. (^|). For g r < and gi = solid lines indicates stable 
radii and dashed line unstable ones. Inset: The eigenvalue A 
vs. g for real VGLE. 

The bifurcation at q = q c is subcritical for g r > 
but supercritical for g r < 0. In the latter case, it gener- 
ates a stable solution with finite R. Under these condi- 
tions, the dependence of the radial velocity v r on R has 
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a zero at some equilibrium distance, as shown in Figs. 3, 
4. For g r > this bound state is unstable. As q fur- 
ther increases, a new pair of stable/unstable bound states 
emerges as a result of a saddle-node bifurcation. This is 
shown in Fig. 4. The number of equilibrium radii is fi- 
nite since the asymptotic behavior of the radial velocity 
is always monotonic. Certainly, only the bound state 
with the smallest radius is important, since at large R 
the "binding strength" decreases. 

The analytical results were compared with numerical 
simulations of Eq. <^). We used the Crank-Nicholson 
method in the domain of 250 x 250 units with 500 x 500 
grid points and with non-reflecting boundary conditions 
(see for details |b| ) . The analytical and numerical results 
appears to be in excellent agreement, as seen in Fig. 3. 

Defect pair in the same field. The interaction of 
oppositely-charged defects in the same field can be re- 
duced to interaction of a single defect with a plane bound- 
ary ||. The problem of interaction in the scalar CGLE 
is determined by the growing exponential solution of the 
stationary linearized problem (see for details Ref . || ) . As 
a result, the velocity due to interaction with the bound- 
ary is given by v ~ e~ pR , where R is the distance to the 
boundary and p is the root with the positive real part 
of the corresponding characteristic equation obtained at 
r 1. As shown in Ref. ||, p is real for q < 0.85 
and complex for q > 0.85. Thus, interaction of the defect 
with a plane boundary (or two defects) is oscillatory, and 
a variety of (unstable) bound states is possible. 



we expect exponentially decaying interaction of the de- 
fects in CVGLE, v ~ e~ pR . For g — ► the corresponding 
exponents should be close to those of the scalar CGLE. 
This exponential interaction is not surprising, since the 
waves emitted by the defects collide and form shocks 
which screen the cores of the defects. 

This exponential interaction is verified by numerical 
simulation of CVGLE. As one sees in Fig. 5, the depen- 
dence of the radial velocity (i.e. the component along 
the line connecting the cores) vs. R is very well approxi- 
mated by an exponential dependence with the exponent 
derived from the solution of the stationary problem. 

We have shown that the interaction between two well- 
separated defects in different complex field is always long- 
range. In a certain parametric domain, these defects may 
form stable rotating bound states. In contrast, the inter- 
action between defects in the same field falls off expo- 
nentially. Our analytical results are limited to a small 
coupling constant g. Additional core instabilities may be 
encountered as g grows M. 
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FIG. 5. The radial velocity v r vs. distance to the bound- 
ary R for q = 1.5 and g = 0.04. The dashed line with circles 
shows the results of simulations with Eq. (^), the solid line is 
a fit by v ~ e~ pR + c.c, where p = 0.746 + i0.94 is the root of 
the characteristic equation of the stationary problem M. 

The interaction problem in the framework of CVGLE is 
similar to that of the scalar CGLE. Looking at station- 
ary perturbations of CVGLE w±, we find that generic 
solutions grow exponentially away from the defect core 
w± ~ Aoe pr , where p is the corresponding root of the 
characteristic equation and Aq is the eigenvector. Thus, 



[1] M.C. Cross and P.C. Hohenberg, Rev. Mod. Phys. 65, 
851 (1993). 

[2] L.M. Pismen, Vortices in Nonlinear Fields, Oxford Uni- 
versity Press (1999). 

[3] L. Gil, Phys. Rev. Lett. 70, 162 (1993). 

[4] L.M. Pismen, Phys. Rev. Lett. 72, 228 (1995); Physica 
D 73, 244 (1994). 

[5] M. Haelterman and A. P. Sheppard, Phys. Rev. E 49, 
4512 (1994). 

[6] T.J. Bridges, Physica D 57, 375 (1992). 

[7] M. San Miguel, Phys. Rev. Lett. 75, 425 (1995). 

[8] E. Hernandez-Garcia, A. Amengual, R. Montagne, M. 
San Miguel, P. Colet, and M. Hoyuelos, Europhys. News 
29, 184 (1998). 

[9] I.S. Aranson, L. Kramer and A. Weber, Phys. Rev. E 47, 
3221 (1993); ibid 47, 4337 (1993) 
[10] P.S. Hagan, SIAM J. Appl. Math. 42, 762 (1982). 
[11] The value of b should not be close to the limit of spi- 
ral core instability, see I.S. Aranson, L. Kramer and A. 
Weber, Phys. Rev. Lett. 72, 2316 (1994); 
[12] I.S. Aranson, A.R. Bishop, and L. Kramer, Phys. Rev. E 

57, 5726 (1998). 
[13] H. Chate and P. Manneville, Physica A 224 348 (1996). 



4 



